Electrostatic theory for imaging experiments on local charges in quantum Hall systems 
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We use a simple electrostatic treatment to model recent experiments on quantum Hall systems, 
in which charging of localised states by addition of integer or fractionally-charged quasiparticles 
is observed. Treating the localised state as a compressible quantum dot or antidot embedded in 
an incompressible background, we calculate the electrostatic potential in its vicinity as a function 
of its charge, and the chemical potential values at which its charge changes. The results offer a 
quantitative framework for analysis of the observations. 

PACS numbers: 73.43.Cd, 73.21.La, 73.23.Hk 
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INTRODUCTION 



Recent imaging experiments 1,2 on quantum Hall sys- 
tems have resolved individual localised states and identi- 
fied discrete charging events in which the charge of these 
states jumps when the mean electron density of the sys- 
tem is altered. The experiments use a scanning probe 
carrying a single-electron transistoBBi to measure locally 
the electrostatic potential and the compressibility of 
the two-dimensional electron gas forming the quantum 
Hall system. For a system close to an integer quan- 
tum Hall plateau, the charging events are believed to 
involve the addition or removal of a single electron 1 , to 
or from the localised state. Close to a fractional quan- 
tum Hall plateau, by contrast, the observed jumps in lo- 
calised charge correspond to the movement of fractionally 
charged quasiparticles^ The latter measurements there- 
fore provide a very direct probe of these quasiparticles, 
whose existence is central to the theory of the fractional 
quantum Hall effect £ 

In this paper, motivated by these experiments, we set 
out a simple description of a localised state in a quan- 
tum Hall system for a regime where behaviour is dom- 
inated by Coulomb interactions. We treat interactions 
using the Thomas-Fermi approximation, making use of 
the well-established picture for screening in integer quan- 
tum Hall systems; 7 : 8 : 9 : 10 : 11 : 12 ! 13 : 14 : 15 in which the sam- 
ple is divided into compressible regions where the local 
Landau-lcvcl filling factor is non-integer, and incompress- 
ible regions where the filling factor is integer. Taking this 
approach, a localised variation in charge density, embed- 
ded in an incompressible background, may be induced 
around a maximum or minimum in the electrostatic po- 
tential due to donors and impurities. In this way, a quan- 
tum dot or antidot is formed with a net charge that is 
an integer multiple of the electron charge for an incom- 
pressible background with integer filling factor. To treat 
localised states in fractional quantum Hall systems, we 
simply assume that quasiparticle charge replaces electron 
charge. While the theory of such quantum dots has been 
discussed in some detail previously; 9 ! 13 : 14 and reduces to 
a standard problem in electrostatics^ a calculation of 



quantities relevant for imaging experiments has not, so 
far as we know, been presented previously. We hope that 
the results we describe here will be useful in further anal- 
ysis of the observations. 



II. MODELLING 

To be definite, we discuss electrons partially filling a 
Landau level to form a quantum dot, which has charge 
density er(r) as a function of position r in the plane of 
the two-dimensional electron gas. An impurity potential 
Vi m p(r) and the screened potential V scr (r) are related by 
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(1) 



where we denote the electron charge by — e. Using the 
Thomas-Fermi approximation for a quantum Hall sys- 
tem, <j(r) =0 in the incompressible region surrounding 
the dot. Throughout the compressible region that makes 
up the dot, screening is perfect and electrons are free to 
adjust their density so that V scr (r) = fi, the chemical 
potential. The screening charge density is restricted to 
lie within the limits < — cr(r) < <y max , where u max is 
the magnitude of the charge density in a filled Landau 
level; we consider only impurity potentials flat enough 
that this upper limit can be ignored. We choose Vi mp (r) 
to have an axially symmetric, parabolic minimum at the 
origin, so that 



V imp (r) =Kr 2 



(2) 



within the radius r ( i of the compressible region. 

Imaging experiments 1 ^* 5 " probe the electrostatic poten- 
tial $(r, z) due to the charge in the localised state rep- 
resented by the dot. We idealise the electron g 
charge sheet of vanishing thickness, located exactly at 
the interface between semiconductor and vacuum, with 
relative dielectric constants S\ and £2 = 1- The resulting 
electrostatic problem is equivalent to one in which there 
is a single medium with dielectric constant e = {e\ + £2). 
This approximation is good provided r ( i is large compared 
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with the thickness of the electron gas and compared with 
its depth below the semiconductor surface, which seems 
to be the case in the experiments of Ref. Q] and 0- The 
potential satisfies Laplace's equation in three dimensions, 
except on the plane of the electron gas within the com- 
pressible region, where the boundary condition 



-e$(r, 0) = (J, - Vi m p(r) for r < r d (3) 

applies. In addition, in the incompressible region, con- 
sistency requires 



-e$(r, 0) > n - V lmp {r) for r > r d . (4) 
The solution can be written in the form 16 



/■oo 

$(r,z)= dkA(k)J (kr)e- k ^, 
Jo 



with 



A(k) 



f(t) cos(kt)dt. 



where, for the parabolic potential of Eq. J2J) , 
2 
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-{2KF - n). 



(5) 



(6) 



(7) 



The charge density in the compressible region is deter- 
mined from 
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dz 



z=0+ 



dz 



(8) 



for r < rd- The value of the chemical potential is fixed by 
the requirements that Eq. is satisfied and that there 
is no divergence in in the charge density: it is 



H = 2Kri 

With this, the charge density for r ^ rd is 



o-(r) 



IQKeeq 



r 2_ r 2 



(9) 



(10) 



(a result given previously in, for example, Ref. IT^) and 
the total charge on the dot is 



Q = 2n / r dr a(r) 



where we have introduced the length scale 
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(11) 



(12) 



It is also useful to calculate the total energy E(Q) of 
the charge on the dot, which can be done by integrating 



the relation \i = —edE(Q)/dQ, using Eqns. © and (fTTf l . 
We find 



E{Q) = ^Ka 2 

5 \ e 



5/3 



(13) 



At this stage, we take account of the fact that charge 
is discrete by setting Q — —Ne, where the number of 
electrons contained in the dot is N = 1, 2, 3, ... In conse- 
quence, the dot radius takes the values 



r d = aN 1/3 . 



(14) 



Having restricted the charge to these discrete values, 
$(r, 0) is no longer related to the chemical potential for 
the sample by Eq. J2J): instead, combining Eq. © and 
Eq. one has 

-e$(r, 0) = 2Krj - V imp (r) for r < r d . 

The values of \x at which charge jumps occur can be found 
by minimising the free energy F — E(Q) — /iN of the dot 
in equilibrium with a charge reservoir, over integer N, 
and considering the result as a function of \x. From the 
expression 



F = -Ka 2 N 5/3 - uN 
5 



(15) 



we find that the values of y, at which the occupation of 
the dot changes between N and N + 1 are 

^ N+1 = b -Ka 2 [{N + I) 5 / 3 - N 5 / 3 ] . (16) 
5 

Next we evaluate the electrostatic potential $(r, z). 
Combining Eqns. ©, Q and iflljl. we have 



AK 



$jv(r,z) = -— / dt(a 2 N 2/3 - t 2 ) x 
J 

x / e~ fe|z| cos(kt)J (kr)dk . (17) 



In this expression, the integral on k can be evaluated 
analytically but the one on t must be done numerically. 
The result can be written in terms of the scaled variables 
p = r/a, ( = z/a and t = t/a as 



with 



<!>N(ap, a() = 



f n (p,0 



Fn(p,Q = ^ dT{N 2 ' 3 ~T 2 ) 



(18) 



1 1/2 



y/X* + 4t 2 C 2 + A 2 
2A 4 + 8r 2 C 2 ) 



where A 2 = ( 2 + p 2 — r 2 . 

Far from the dot, for (p 2 + C 2 ) >> iV 2/3 , these expres- 
sions reduce to 



$Ar(r,z) 



-Ne 1 



47T££ yjr 2 + Z 2 



(19) 
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FIG. 1: The scaled electrostatic potential Fisr(p, f) as a func- 
tion of scaled radius p and height £ from the centre of the dot, 
for N = 1 (lower surface), N = 2 and iV = 3 (upper surface). 
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FIG. 2: _Fjv as a function of p, at fixed height £ above the 
dot. In each graph, the curves are for N = 1 (lowest), N — 2, 
N = 3, and AT = 4 (highest). 



as expected. 

The dependence of the function Fjv (p, C) 011 P an d C is 
illustrated in Fig. ^ and its variation with p at fixed £ is 
shown in Fig. [21 

Since in experiment this potential will add to other 
contributions, for example, from fixed background 
charges, it is useful to focus on the potential changes aris- 
ing from jumps in the charge of the dot. These changes 
are proportional to AF N (p,() = F N+1 (p,() ~ F N (p,Q, 



and this function is shown in Fig. 

AF N (p,£ = 0) 
1 




AF N (p,£ = l) 
0.8 




FIG. 3: Electrostatic potential changes (proportional to 
AFjv) produced by jumps in the charge of the dot, as a func- 
tion of p, at fixed distances, £ = and £ = 1 above the dot. 
In each graph, the different curves are for N — 1 (the highest 
curve at p = 0), N = 2, N = 3, N = 10 and N = 20 (the 
lowest curve at p = 0) . 



III. DISCUSSION 

There is scope to compare these results with experi- 
ment in several different ways. 

First, the most striking feature of the observations is 
the fact that, considering behaviour as a function of av- 
erage electron density n and flux density B, a particular 
charging event takes place on a line in the n-B plane 
which is parallel to one of the lines of integer filling fac- 
tor v. Such behaviour is built into the model we have 
studied. In particular, suppose that Nl Landau levels in 
the sample are completely filled, so that the charge den- 
sity within the quantum dot (c(r) in Eq. Q) lies in the 
{Nl + 1) th level. In that case, if n and B vary together 
along a line in the n-B plane parallel to v — Nl, the 
charge density variation in the sample is uniform in space 
and the screened potential remains constant; as a con- 
sequence, the charge of the quantum dot is unchanged. 
Conversely, charging events are produced by moving in a 
perpendicular direction in the n-B plane. This account 
omits the single-particle contribution, (Nl + l/2)hw c , to 
the energy of the charge within the dot. The approxi- 
mation is justified because the electrostatic part of the 
energy of the two-dimensional electron gas is dominant. 



4 



Turning to more specific aspects of our modelling, it 
is useful to focus on results that are independent of the 
model parameter K in Eq. (j2J and of measurement cal- 
ibration. Two such results (which are physically related 
to each other) are the power laws appearing in the de- 
pendence of dot radius on electron number, Eq. 1|14|) . and 
in the chemical potential values at which charge jumps 
occur, Eq. (jlbf) . While dot radius is probably difficult 
to measure precisely, because of issues of resolution, as 
discussed below, the relative size of chemical potential 
steps required to add a sequence of charges should be an 
accessible quantity. Specifically, charge jumps are pro- 
duced experimentally by a change in the backgate voltage 
applied to a sampleti*^ fitting the ratio of voltage steps 
for successive jumps to Eq. Qlfijl would provide a test of 
the theory we have presented and a determination of the 
number of electrons within the dot. Deviations from the 
theory would arise either if the confining potential is not 
parabolic, or, more interestingly, if many-body correla- 
tions within the compressible region, which are omitted 
from Thomas-Fermi theory, make an important contri- 
bution to the total energy of the electrons in the dot. 
Even in these cases, we expect as a robust feature a de- 
crease in the size of voltage steps between charge jumps 
as electron number increases. 



In addition, one can attempt an absolute comparison of 
theoretical and experimental quantities. As an illustra- 
tion, suppose a — 200 nm and N = 10, so that Td = 430 
nm, and consider a measurement of the potential by a 
scanning probe at a height z — 200 nm above the sam- 
ple. Taking, for GaAs, e\ — 13, we find from Eq. (|18J) a 
change in electrostatic potential when a further electron 
is added, of size A$i (r = 0,z= 200nm) = 520ajV. This 
is similar to the step size of 180/xV reported in Ref.H an 
exact match could presumably be arranged by adjusting 
a, N or z. Beyond this, one can regard our calculation 
of AF(p,Q, illustrated in Fig. |3 as a determination of 
the resolution function for the imaging technique. 

In summary, we have presented a simple model for the 
imaging experiments of Ref . Q] and 0- A closer compar- 
ison between observations and calculations should help 
determine the numbers of electrons contained in localised 
states and the spatial size of these states, while deviations 
of measurements from this theory may be an indication 
of correlation effects. 

We are grateful to Amir Yacoby for discussions and 
for preprints of Refs. Q] and U The work was sup- 
ported in part by CAPES, and by EPSRC under Grant 
GR/R83712/01. 
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